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We interpret the unique high-energy spectrum of the millisecond pulsar TfSR J0218+4233i 
within polar cap scenarios. We show that the spectral data from BeppoSAXliil and EGRET 
impose very restrictive limitations on possible radiation mechanisms, energy spectrum of radi- 
ating charges as well as viewing geometry. Theoretical spectra are able to reproduce the data, 
however, this can be achieved provided very special - unusual within the conventional polar 
cap picture - conditions are satisfied. Those include off-beam viewing geometry along with 
one of the following alternatives: f ) strong acceleration of secondary pairs; 2) broad energy 
distribution of primary electrons extending down to 10 5 MeV; 3) high-altitude synchrotron 



1 Introduction 

PSR J0218+4232, with the spin period P = 2.3 ms and the inferred dipolar magnetic field at 
polar cap B pc ^_8-6 x 10 8 G, is the only millisecond pulsar which has been marginally detected 
above 100 MeVEa. A broad-band high-energy spectrum of this object exhibits unusual features, 
completely different from those observed among young gamma-ray pulsars: Above 100 MeV the 
photon index a p h ~ —2.60 and the spectrum resembles the very soft spectra of middle latitude 
unidentified EGRET sourcesi Within the BeppoSAX range the spectrum is extremely hard: 
a ph ~ -0.61 ± 320 

Kuiper et alp noticed that neither the polar cap nor the outer gap models could naturally 
explain the spectrum. The aim of this work is to carefully examine an ability of polar cap models 
to reproduce the spectrum of J0218+4232. 



2 Directional characteristics of pulsar spectra predicted by the polar cap model 

As we show in the acompanying paper (Wozna et al., these proceedings) the viewing geometry 
effects have crucial significance for the appearance of a pulsar spectrum for a given observer. For 
purely dipolar magnetic field the geometry is determined by two angles: a - the angle of dipole 
inclination relative to the rotation axis, and £ - the angle between an observer's line of sight and 
the rotation axis. Various combinations of a and £ result in a large variety of spectral shapes 
and pulse profiles in high-energy domain. However, a particular energetic history of electrons 
in pulsar magnetosphere enables to extract two main cases: the on-beam and the off-beam 
geometry. 

Fig. la presents the Lorentz factor 7 of electrons injected at the surface of neutron star along 
the "last open" magnetic field lines as a function of angle 6 m between magnetic dipole axis and a 



Figure 1: a) Electron energy 7 as a function of magnetic co- 
latitude of photon propagation direction # m for P = 23 ms and 
70 = 2 x 10 7 ; b) Characteristic energy of curvature photons e cr 
as a function of # m for the electron energy shown in the upper 
panel. The dashed region (SR) indicates the range of 6 m where 
most of synchrotron photons are emitted. The figure implies that 
for most viewing angles £ one of two main shapes of spectrum can be 
recorded: In the case of the on-beam geometry the line of sight sam- 
ples 9m < 1.5 0p C + 1/711 — 32° and the observed spectrum consists of 
curvature and synchrotron components. The first one extends up to 
~ 10 5 MeV with a p h ~ —1.66 and the latter one, with a p h = —1.5, 
dominates below a few tens of MeV. For the off-beam geometry only 
the curvature component can be detected with a p h = —2/3 below 
~ 10 MeV and a p h > 2 within the range of its high-energy cutoff 
around 100 MeV. PSR B1821-24 and PSR J0218+4232 seem to be 
the examples of these two cases, respectively. An intermediate case 
- with the line of sight just grazing the heighest-energy gamma-ray 
beam - is less probable. The gray band marks the range of 8m sam- 
pled by the line of sight at £ = 48° for a — 8°. The spectrum 
corresponding to this off-beam viewing geometry is shown as a thick 
solid line in Fig. 2. 

local tangent to magnetic field line at the electrons' position. The value of 6 m = 26. °4 ~ 1.5 9 pc 
corresponds to the direction of B at the polar cap rim (9 pc is the angular radius of the polar 
cap, measured from the center of neutron star). The rotation period P = 2.3 ms and the initial 
Lorentz factor 70 = 2 x 10 7 were assumed in these calculations. The electron energy losses 
noticeable in Fig. la are purely due to the emission of curvature radiation]] The CR energy loss 
rate of an electron is initially huge, but because of its strong dependence on electron energy 
(7cr oc 7 4 //jp r , p cr is the local radius of curvature of magnetic field lines) it becomes negligible 
after the electron traverses a length which is small in comparison with the light cylinder radius 
R\ c . Thus, after the initial rapid drop, the electron's energy starts to decrease very slowly 
(Fig. la). 

The characteristic energy 6 cr of curvature photons emitted by the electron in different direc- 
tions 6 m is shown in Fig. lb. Because e cr is very sensitive to 7 (e cr oc 7 s / p cr ) the high-energy cut- 
off in the CR spectrum decreases rapidly for increasing angles 9 m but then it starts to approach 
"asymptotically" the value of ~ 30 MeV. This is in part due to the slower decrease in 7 at higher 
altitudes and rn part due to an increase in p cr which starts to take place for # m > 68°. In fact, 
Rudak & DyksEJ (RD99) estimated an "absolute minimum" of 7t, re ak — 3.8x 10 6 (p cr /(10 7 cm)) 1 / 3 
for the Lorentz factor of CR-cooled electrons (eq. (7) in RD99) and the absolute lower limit of 
~ 150 MeV for e cr (eq. (8) in RD99). These values are in reasonable agreement with the exact 
results presented in Fig. 1, given the crude method used in RD99. We recall here that the lower 
limit for e cr practically does not depend on any pulsar parameters (eg. exact calculations for 
rotation period P = 0.1 s g 1V6 €cr ^ 20 MeV at R Xc ). 

An important implication of Fig. lb is that pulsar radiation pattern can be considered as 
consisting of two components: a hollow cone of very high-energy emission extending up to GeV 
range (with opening half- angle 8 m ~ 1.5 9 pc ) and a much less anisotropic emission peaking 
close to 30 MeV. When the line of sight crosses the hollow cone beam (on-beam geometry) 
the recorded spectrum consists of two components: a curvature component and a synchrotron 
component due to synchrotron radiation (SR) from secondary electron-positron pairs. The 
shape of this spectrum is similar to the total spectrum emitted in all directions, shown in fig. 1 
of RD99. Its qualitative features can also be assesed from Fig. lb: The spectrum extends 




"Resonant inverse Compton scatterings cannot influence the election's energy due to small relative energy 
losses per scattering in a weak magnetic field (eg. Dyks, Rudak, Buliklj) whereas non-resonant scatterings (they 
occur in the Klein-Nishina regime) are too rare for surface temperature values relevant for neutron starsu 
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Figure 2: Comparison of theo- 
retical pulsar spectrum for the 
off-beam geometry (thick solid 
line) with the spectrum observed 
for PSR J0218+4232 (diamonds 
- BeppoSAX, triangles - OSSE, 
squares - COMPTEL, circles - 
EGRET). The BeppoSAX points 
suggest a presence of additional 
spectral component in the X-ray 
domain. Such a component would 
arise due to CR of secondary e 
pairs, provided that at least a few 
hundred of them per primary are 
first accelerated up to 7 ~ 10 s 
MeV (see Section 3). 



up to GeV range and is relatively soft within the EGRET range because it is composed of 
instantaneous CR spectra with different values of e cr (Fig. lb). The value of photon index a p h 
in this spectral range depends on the viewing angle £ as well as on the rotation period P. For 
millisecond periods it is close to a p h ~ —5/3. Within the entire X-ray range (0.1 keV - a few 
MeV) the on-beam spectrum is dominated by the synchrotron component with the well-known 
photon index a p h = —1.5 (RD99). 

As noted by Kuiper et ali, these on-beam characteristics are in clear disagreement with 
the high-energy data on J0218+4232. We, therefore, propose the off-beam geometry for this 
pulsar ie. the case when the line of sight misses the hollow cone beam. With this assumption 
the problem of wrong slopes within the BeppoSAX and the EGRET range can be solved. 
Since CR photons of the highest-energy are converted into the e ± pairs with Lorentz factors 
7|| ~ 10 (in B ~ 10 9 G), most of synchrotron photons follow the direction of the "parent" 
CR photons and are constrained to a narrow range of angles around 1.5 9 pc with a spread of 
1/7|| ~ 0.1 sr (Fig. lb). Therefore, in the off-beam case one misses the SR emission and the 
spectrum consists of the CR component only. Thick solid line in Fig. 2 presents this kind of 
spectrum calculated for a = 8° and £ = 48° (which corresponds to 40° < 9 m < 56° in 
Fig. 1). The spectrum resembles very closely the well-known instantaneous CR spectrum due 
to monoenergetic electrons with no cooling. The reason for this is clearly shown in Fig. lb: 
because of the near alignment of the dipole and the rotation axis (a ~ 8°) the off-beam line of 
sight samples only a very limited range of 9 m : between 40° and 56° (grey band in Fig. lb). U 
Below a high-energy cutoff at ~ 100 MeV the CR spectrum has the photon index a p h = —2/3. 

If the X-ray emission from J0218+4232 is generated close to the neutron star surface, it can 
be explained only as a pure CR emission (note that a p h = — 2/3 matches well the BeppoSAX 
photon index —0.61 ±0.32). Any SR emitted close to the surface would have a broken power-law 
shape with a p h = —1.5 within the X-ray range, in clear disagreement with the BeppoSAX data. 
The slope a p h = —2/3 of the CR component does not depend on the viewing angle £ as long as 
the line of sight misses the SR component, (which is concentrated around 1.5 9 pc ), ie. as long as 
C > a + 1.5 6 pc + 1/7|| ~ 40° (see Wozna et al., this volume). 

Furthermore, owing to the off-beam geometry the high energy cutoff in the total spectrum 
(which now consists of just a CR component) occurs at a relatively low photon energy of ~ 100 
MeV. The cutoff is not caused by magnetic absorption - it corresponds to a maximum energy of 

b Interestingly, the shape of the spectrum would not change much even if the line of sight sampled much larger 
range of # m within the off-beam region (say between 40° and 140°). This is because of the remarkable stability 
of e cr at high altitudes - see Fig. lb. 



those electrons which emit observable CR. This is why the cutoff's shape can easily mimic the 
very soft spectral shape as suggested by the two EGRET points (Fig. 2). Harding & Zhang§ 
used similar viewing geometry arguments to explain the very soft spectra of unidentified EGRET 
sources. 

However, with the off-beam CR spectrum normalized to reproduce the level of the EGRET 
emission, the BeppoSAX data points are located about 3 orders of magnitude above the ex- 
trapolated level of the CR component (thick line in Fig. 2). The intrinsically hard shape of the 
off-beam spectrum results in too low ratio of X-ray level to gamma-ray level of emission. In our 
opinion it is not possible to solve this problem within the standard polar cap model. 

3 Can the secondary pairs emit CR photons within the keV range? 

In order to explain the BeppoSAX data we propose a contribution of CR from secondary 
pairs accelerated within the polar gap. Upon the acceleration energy 711 of the pairs should 
reach values in the range between 1.5 x 10 5 and 6 x 10 5 . For 711 < 7b re ak ~ 3 x 10 6 the CR energy 
loss length scale considerably exceedes the size of magnetosphere (R\ c )- The pairs do not suffer 
thus considerable energy losses and at all altitudes emit roughly the same spectrum of CR. In 
Fig. 2 it is shown as a thin solid line (we took 711 = 3 x 10 5 ). 

The level of CR spectrum at energies below the characteristic energy e cr oc 7 3 //9 cr does not 
depend on 711, but solely on the radius of curvature of magnetic field lines p cr . Since the pairs 
follow nearly the same field lines as the primary electrons, the ratio betwen the level of BeppoSAX 
data points (diamonds in Fig. 2) and the level of X-ray emision from primary electrons (thick 
solid line) is a direct measure of the number of pairs n± per single primary electron required to 
reproduce the level of the X-ray data. The level of the primary CR spectrum (thick solid line) 
within the BeppoSAX range can be determined by fitting the high-energy cutoff of the spectrum 
to the EGRET data points. Apart from n±, the same fit could also determine the value of £, 
because the shape of the cutoff depends on the viewing angle .fl, however, the EGRET data do 
not allow to determine n± and £ with high accuracy. "By eye" fits give n± of at least a few 
hundred; for ( = 48° (Fig. 2), n± ~ 10 3 . 

Thus, to reproduce the broad-band high-energy spectrum of J0218-I-4232, at least a few hun- 
dred of pairs should acquire an energy roughly equal to 1% of primary electron energy. According 
to most works on the physics of the pair formation front (eg. Harding & MuslimovB) acceleration 
of such a large number of pairs within the polar gap is difficult, because redistribution of the 
pairs screens out the electric field of the polar gap. However, existing works on the physics of 
the polar gap neglect the fact that the pairs are created at different magnetic field lines than 
those along which parent primary electrons propagate. Because of this ID-approach as well as 
plenty of other approximations our present knowledge of processes at the pair formation front 
is far from being well established. 

Because of the bimodal energy distribution of radiating charges (primary electrons, secondary 
pairs), the CR spectral components which correspond to them are separated by a dip. A precise 
shape of this dip would depend on the shape of high-energy cutoff in the energy distribution of 
e pairs and the dip may extend in principle between ~ 10 keV and 100 MeV (dashed line in 
Fig. 2 suggests a possible situation). 

4 What is the energy spread among primary electrons? 

One of popular assumptions present in polar cap models is that all primary electrons accelerated 
within the gap acquire roughly the same energy. This is not realistic because the potential drop 

c Because of the near alignment of the magnetic dipole, the value of £ cannot be well determined from radio 
polarisation data on J0218+4232 (Stairs et altJ) 



across the Dolar cap depends on magnetic colatitude: it decreases towards the rim of the polar 
cap (eg. I'll). Moreover, stochastic energy losses caused by resonant Compton scatterings may 
significantly soften the energy spectrum of primary electrons provided the surface temperature 
T and the magnetic field B are sufficiently high (see fig. 4 in paper i). In the case of weak-1? 
millisecond pulsars the latter mechanism is not efficient because the energy loss rate due to the 
resonant ICS (oc B) is low. However, no calculations of primary electron energy distribution 
have been performed so far for very high temperatures in excess of 10 7 K - a value considered 
for JO218+42320 

Suppose that the primary electrons assume the steady-state energy distribution in a form 
of a power law dN e /d'y oc 7 P extending down to 7 m i n = a few x 10 5 . Then, the resulting CR 
spectrum could easily reproduce both the observed slopes and levels of X-ray and gamma-ray 
emission for appropriate values of p and 7 m i n . It would have a broken power law shape with 
a p h = —2/3 within the BeppoSAX range, and a p h = (p — l)/3 above a break somewhere 
between 10 keV and 100 MeV (the exact value of the break energy et, r = e cr (7 m i n )). Again, if 
primary electrons really achieved 7 in excess of 10 7 in the polar gap of millisecond pulsars, the 
agreement of the model spectrum with EGRET data points could be easily achieved only by 
a proper choice of the viewing angle £ for the off-beam geometry. There would be no dip in 
the spectrum and a measure of a p h and et> r would give us direct information about p and 7 m i n , 
respectively. BeppoSAX data points and COMPTEL upper limits constrain a p h to the range 
between —1.5 and —2.15 which corresponds to rather soft electron energy distribution with the 
index —3.5 > p > —5.5 and with the low energy cutoff at 10 5 < 7 m i n < 2.5 x 10 5 . 

5 High-altitude synchrotron emission 

In the case when the synchrotron radiation (SR) is emitted close to the neutron star surface, its 
spectrum extends down to a blue-shifted local cyclotron energy j\\tiujB with a p h = —1.5 (RD99). 
For high-altitude emission, however, the cooling length scale due to SR can become longer than 
the length scale of decrease in B. Simple calculations show that this may happen for a local 
magnetic field B < 2 x 10 6 G. For PSR J0218+4232 such a field is expected at r > 0.5 R ic . 

In the case of SR in such a low B, electrons do not loose their entire energy ^y±mc 2 corre- 
sponding to their motion across B. In consequence, a new break appears in the SR spectrum, at 
photon energy greater than 'y^huiB (Chang et alp). Below the break, the photon index assumes 
the value a p h — —2/3, characteristic for instantaneous spectrum of SR (or CR) emission. 

Fig. 3 shows that such a kind of spectrum is also able to reproduce the high-energy data on 
J0218+4232. Between the BeppoSAX and EGRET range the SR spectrum has a well known 
slope of a p h = —1.5, and its flux level is about 10 times below the upper limits from OSSE 
and COMPTEL. This SR spectrum was calculated within the following simplified model: The 
electrons were injected at a distance of 0.8i?i c from a neutron star, with initial energy 7 = 3 x 10 7 
and were propagated radially up to f>R\ c with a constant value of 711 = 10 5 (ie. constant pitch 
angle during the SR emission was asumed). Viewing geometry effects were not taken into account 
in these simple calculations but the lack of strong GeV emission from J0218+4232 again implies 
the off-beam case. The agreement between the data and the model spectrum of SR shown in 
Fig. 3 is remarkable, however, the emission of 100 MeV synchrotron photons close to the light 
cylinder requires very high energy of primary electrons near R\ c : (7j_ > 10 2 , 711 ~ 10 5 ). Various 
mechanisms of acquiring 7j_ > 1 in the outer parts of magnetosphere have been considered so far 
within the polar cap model (eg. Malov & Machabeli0) but they usually concerned the optical 
emission from pulsars. 




Figure 3: Comparison of theoret- 
ical pulsar spectrum for a pure 
SR emission close to (and beyond) 
the light cylinder (thick solid line) 
with the spectrum observed for 
PSR J0218+4232. Note that the 
characteristic slope of SR emis- 
sion with a p h = —1.5 for e > 
30 keV changes at 30 keV into 
a pb ~ —2/3 within the Bep- 



poSAX range. 
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6 Conclusions 

The combined X-ray and gamma-ray data on the millisecond pulsar J02 18+4232 impose severe 
limitations on possible radiation mechanisms, energy spectrum of radiating charges as well 
as viewing geometry. The standard polar cap model is not able to explain the spectrum of 
J0218+4232 even when viewing geometry effects are taken into account. We find three possible 
interpretations which require non-orthodox assumptions about the electron energy distribution 
or emission altitude. The spectra corresponding to these possiblities have unique features (the 
MeV dip or the characteristic slope of SR) which may enable to identify them with a high- 
sensitivity and high-angular resolution gamma-ray telescope (INTEGRAL?). 
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